Half-metallic properties of atomic chains of carbon-transition metal compounds 
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We found that magnetic ground state of one-dimensional atomic chains of carbon-transition metal 
compounds exhibit half-metallic properties. They are semiconductors for one spin-direction, but 
show metallic properties for the opposite direction. The spins are fully polarized at the Fermi 
level and net magnetic moment per unit cell is an integer multiple of Bohr magneton. The spin- 
dependent electronic structure can be engineered by changing the number of carbon and type of 
transition metal atoms. These chains, which are stable even at high temperature and some of which 
keep their spin-dependent electronic properties even under moderate axial strain, hold the promise 
of potential applications in nanospintronics. 

PACS numbers: 61.46.+W, 72.25.-b, 75.50.Cc, 71.20.Be 



Spin-dependent electronic transport has promised rev- 
olutionary applications using giant-magneto-resistance in 
magnetic recordingand nonvolatile memories Q, 0, Q- 
Half-metals (HM) ]1 are a class of materials, which 
exhibit spin-dependent electronic properties relevant to 
spintronics. In HMs, due to broken spin-degeneracy, en- 
ergy bands E n (k, f) and E n (k, j) split and each band ac- 
commodates one electron per k-point. Furthermore, they 
are semiconductor for one spin direction, but show metal- 
lic properties for the opposite spin direction. Accord- 
ingly, the difference between the number of electrons of 
different spin-orientations in the unit cell, N — — Ni, 
must be an integer and hence the spin-polarization at 
the Fermi level, P = D{E F ,]) - D{E F ,i)/[D{E F ,} 
) + D(E F , [)], is complete pj. This situation is in con- 
trast with the ferromagnetic metals, where both spin- 
directions contribute to the density of states at E F 
and spin-polarization P becomes less than 100%. Even 
though three-dimensional (3D) ferromagnetic Hcuslcr al- 
loys and transition-metal oxides exhibit HM properties 
[a] , they are not yet appropriate for spintronics because of 
difficulties in controlling stoichiometry and the defect lev- 
els destroying the coherent spin-transport. Zinc-blende 
(ZB) HMs with high magnetic moment fi and high Curie 
temperature T c > A00K (such as CrAs, and CrSb in ZB 
structure) have been grown only in thin- film forms 0. 
More recently, it has been predicted that four new ZB 
crystals can be HM at or near their respective equilib- 
rium lattice constants @. 

In this letter, we report that very simple and sta- 
ble one-dimensional (ID) structures, such as linear 
atomic chains of carbon-transition metal compounds, i.e. 
C„(TM), show half- metallic properties. The prediction 
of half-metallic behavior in ID atomic chains is new 
and of fundamental interest, in particular in the field of 
fcrmionic excitations with spin degree of freedom. Be- 
sides, the present finding may lead to potential appli- 
cations in the rapidly developing field of nanospintron- 



ics, such as tunnelling magnetoresistance, spin valve, and 
nonvolatile magnetic devices. 

In earlier transport studies, the spin-direction of con- 
duction electrons was generally disregarded, in spite of 
the fact that spin orientation of electrons decays much 
slower than their momentum [3j. Magnetic ground state 
of TM-adsorbed single-wall carbon nanotubes (SWNT) 
lioj l. spontaneous spin-polarized electron transport 
through finite TM wires [llj], and oscillatory spin- 
polarized conductance and spin-valve effect through fi- 
nite carbon wires capped with Co atoms in between gold 
electrodes have been treated recently. However, half- 
metallicity predicted in periodic C„(TM) is a behavior 
fundamentally different from those magnetic properties 
found in earlier systems in Ref. 0, 0, fill and is 
a novel feature of ID systems that allows both semicon- 
ducting and metallic properties coexisting in the same 
structure. 

Our predictions are obtained from the first-principles 
pseudopotential plane wave calculations within DFT us- 
ing generalized gradient approximation (GGA) and ul- 
trasoft pseudopotentials 13]. C n (TM) chains have been 
treated in supercell geometry, and Brillouin zone is sam- 
pled by 10 — 80 special k-points depending on the super- 
cell size. All the atomic positions, and the supercell lat- 
tice parameter c along chain axis are optimized by min- 
imizing the total energy, E^ , the forces on the atoms, 
as well as the stress of the system calculated by spin- 
polarized (or spin-relaxed) GGA. Since AE = E™ — E^ p 
(where EjP is spin-unpolarized total energy) and net 
magnetic moment /x are both positive, these compound 
chains have ferromagnetic ground state [l4|. While our 
study has covered a large family of C„(TM) chains, our 
discussion will focus on C„Cr. 

Finite size linear chain of single C atom has already 
been synthesized experimentally (lifr . The double bond 
between carbon atoms and doubly degenerate 7r-band 
crossing the Fermi level underlie the stability of the chain 
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FIG. 1: Transition state analysis for two different reaction 
paths (a) Path-1: Variation of energy, Et(c), for the transi- 
tion from the Cs+Cr initial state to the linear CsCr HM as 
final state: Qi-tf and Q/^i are energy barriers involved in 
the transitions. Zero of energy is taken relative to the free Cr 
atom and free periodic C-linear chain, (b)-(d) Path-2: Varia- 
tion of interaction energy with distance d between C7 and Cr, 
CrCr and C, CrCrC and C for adatom attaching from one 
end. 



and its unusual electronic properties. We examine the 
formation of linear C5& HM by performing transition 
state analysis along two different reaction paths. Nor- 
mally, a Cr atom is attracted by C-linear chain and even- 
tually forms a bridge bond over a C-C bond. We take 
this bound state (specified as Cs+Cr) as an initial state 
of the first reaction path for transition to the final state 
corresponding to the linear C5& HM as illustrated in 
Fig- E a ) uM ■ The energy barrier necessary to go from 
the initial state to the final state is Qi^/=1.86 eV. Once 
the final state has formed, it is prevented from going back 
to the initial state by a significant barrier of Q/^i=1.04 
eV. However, the energy barrier Qi_ >/ disappears totally 
and hence the process becomes exothermic for the second 
reaction path where HM is grown from one end of the 
chain by attaching first Cr atom then C atoms sequen- 
tially. Each ad-atom (Cr or C) is attracted to chain and 
eventually becomes bound to it as shown in Fig.^b-d). 
This analysis lets us believe that half- metallic C n (TM) 
chains are not only of fundamental interest, but also can 
be realized experimentally. 

Whether a periodic C„Cr linear chain is stable or it 
can transforms to other structures has been examined by 
an extensive investigation of Born-Oppenheimer surface. 



TABLE I: Results of spin-polarized first-principles calcula- 
tions for C„Cr linear chains. AEt is the difference between 
spin-paired (non-magnetic) and spin-polarized (magnetic) to- 
tal energies, c is the optimized ID lattice parameter, fi is the 
total magnetic moment per unit cell in units of Bohr magne- 
ton, fiB- M, S, and HM stand for metal, semiconductor, and 
half-metal, respectively. By convention majority and minor- 
ity spins are denoted by f and J.. The numerals in the last 
column are the band-gap energies in eV. 
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Local minima of the total energy have been searched 
by optimizing the structure starting from transversally 
displaced chain atoms for varying lattice parameters. 
The linear chain structure has been found to be sta- 
ble and energetically favorable relative to zigzag struc- 
tures. The phonon calculations of C„Cr, yielding pos- 
itive phonon frequencies (fixo(fc = 0) =89, 92, 411 
cur 1 : Q LO (k = 0) =421, 1272, 1680 cm- 1 for n =3 and 
n TO (k = 0) =13, 71, 353, 492 cm" 1 ; n LO {k = 0) =489, 
1074, 1944, 2102 cm -1 for n =4) corroborate the above 
analysis of stability. However, for n =9 some of the fre- 
quencies get negative indicating an instability for a large 
n. In addition, we performed high temperature (T=750- 
1000 K) ab-initio molecular dynamics calculations using 
Nose thermostat, where atoms are displaced in random 
directions. All these tests have provided strong evidence 
that the linear chain structures with small n are stable. 
To weaken the constraints to be imposed by supercell ge- 
ometry, calculations have been done by using double su- 
percells including two primitive unit cells of the chains. 
Peierls instability that may cause the splitting of metal- 
lic bands at the Fermi level did not occur in C n Cr linear 
chain structures. Table I summarizes the calculated mag- 
netic and electronic properties of C n Cr linear chain. We 
also examined how an axial strain may affect the half- 
metallic behavior of these chains. HM character of C4Cr 
was robust under e = ±0.05. Small band gap CsCr re- 
mains HM for e = 0.05, but is rendered a ferromagnetic 
metal under e = —0.05. While C3& changes to a semi- 
conductor under e = 0.10, it becomes a ferromagnetic 
metal with fi = 3.1 under e = —0.10. 

Spin-polarized electronic band structures are strongly 
dependent on n. For example, all C n Cr we studied are 
HM except CCr, which is a semiconductor. For even 
n, majority spin bands are metallic, but minority spin 
bands are semiconducting with large band gaps (£^4 ^3 
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eV). This situation, however, is reversed for odd n, where 
majority spin bands become semiconducting with rela- 
tively smaller gaps (E g ^ ~0.5 eV), but minority bands 
are metallic. This even-odd n disparity is closely related 
to bonding patterns in different chains. For example, for 
odd n = 3, respective bond lengths are in A -C-1.28-C- 
1.28-C-1.95-Cr-, and for even n = 4, -C-l.25-C-l.33-C- 
1.25-C-2.1-Cr-. It appears that, while double bonds are 
forming between all atoms for odd n, for even n triple 
and single bonds form alternatingly between C atoms, 
and single bonds occur between C and Cr atoms with 
relatively longer bond lengths [l^. Consequently, the 
overlap between Cr and C orbitals and hence relative en- 
ergy positions of bands vary depending on whether n is 
even or odd. The C n Co linear chains exhibit also similar 
spin-dependent electronic properties. While CCo chain is 
a ferromagnetic metal with fx =0.5 lib, C„Co chains (for 
2< n <6 studied in this work) are HM. They have even- 
odd n disparity; E g is relatively larger for even n but 
decreases as n increases. Interestingly, the odd n C n Co 
chains, where nearest carbon atoms have spins smaller 
but opposite to that on Co, exhibit a ferrimagnetic be- 
havior. 

Half-metallic electronic structure and resulting spin- 
dependent properties of C„Cr linear chains are shown by 
the bands and density of states in Fig. El The odd-even n 
disparity is clearly seen. The double degenerate 7r-band 
(denoted by m [ for n =3 or m | for n — 4) is half-filled 
and determines the position of Ep. The band gap of 
semiconducting states, which have spin in the direction 
opposite to that of m-band, occurs between the filled flat 
^i-band and empty conduction ci-band. According to 
these bands, the equilibrium ballistic conductance of the 
infinite CaCr is G [=2e 2 / h for minority spin, but zero for 
majority spin. The calculated spin projected total den- 
sity of states (TDOS) in Fig. shows the energy spec- 
trum of majority and minority spin states in an interval 
±2 eV around Ep. The band gap for one spin direction, 
and finite density of states at Ep for the opposite spin are 
clearly seen. This is a dramatically different finding than 
those of Ref. Orbital projected local densities 

of states at Cr and C atoms show the orbital composition 
of the spin-polarized bands. The m-band is composed of 
Cr-3c? and mainly first-neighbor G-2p orbitals at Ep, that 
is the p — d hybridization. The flat ui-band nearest to 
Ep is derived from the Cr-3d and 4s states. The empty 
ci-band originates from C-2p and Cr-3cZ states. 

The effect of spin-orbit (SO) coupling on the HM prop- 
erties of CaCr has been calculated by using all-electron 
DFT (Wien2K) code. We found the splitting is very 
small and £ T (with SO)-E T (without SO)=-7.9 meV. As 
illustrated in Fig. HJc), the difference between TDOSs 
calculated with and without SO-coupling is negligible 
and hence the effect is not strong enough to destroy 
the half-metallic properties. This conclusion obtained 
for n =3 can apply to other chain structures in Table I, 
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FIG. 2: (a) Energy band structure of CaCr; corresponding 
total density of states (TDOS) for majority (f) and minor- 
ity (|) spins; orbital projected local density of states at Cr 
atom (PDOS/Cr), at C atoms first nearest neighbor to Cr 
(PDOS/C1). (b) C 4 Cr. State densities with s,p,d orbital 
symmetry in PDOSs are shown by thin continuous, dotted, 
broken lines, respectively. Zero of energy is set at Ef- Metal- 
lic band crossing the Fermi level, highest valence and lowest 
conduction bands are labelled by m, vi and ci, respectively, 
(c) Majority and minority TDOS of CsCr calculated with and 
without spin-orbit coupling (SO). 



since the Cr-C interaction decays fast beyond the first 
nearest neighbors of Cr atom. 

Further insight about the electronic structure and 
the character of bonding can be gained by examining 
the charge distributions associated with selected bands 
shown in Fig. |3J The charge density of the metallic spin- 
state at Ep is obtained by averaging charges of states 
of m-band having energy ±0.02 eV around Ep. They 
are formed from the bonding combination (or p — d hy- 
bridization) of 0-2^2;^ and Cr-3d X z,yz valence orbitals of 
CsCr. In the case of even n (n=4) it corresponds to 
an antibonding combination of the above orbitals with 
enhanced Cr-3d contribution. The charge density of v\- 
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FIG. 3: Charge density contour plots of linear chains of CsCr, 
C4& compounds on a plane through the chain axis. p m ,foii is 
the charge density of metallic spin-states within energy range 
i?F±0.02 eV. p vl ,]oTi and p cl ,torj are the charge density of the 
highest valence V\ and the lowest (empty) conduction band 
ci, respectively. pr,i and pr,i are total charge density due 
to majority spin states and minority spin states, respectively. 
Ap is contour spacing. 



band is due to non-bonding Cr-4s — 3d z 2 orbitals in C^Cr. 
For even n case, the C-2p xy contribution is pronounced. 
Charge density of the Ci-band suggests the antibonding 
combination of p — d hybridized states. 

The band structure and charge density plots suggest 
that the p — d hybridization between neighboring C and 
Cr orbitals and resulting exchange splitting of bands in 
different spin directions give rise to the ferromagnetic 
ground state of the above chains. An additional ingre- 
dient, namely cylindrical symmetry of the 7r-bonds in 
the carbon chain provides conditions to result in inte- 
ger number of excess spin in the unit cell which is re- 
quired to achieve half-metallic behavior. The ferromag- 
netic ground state with /1 values integer multiple of [Ib 
per unit cell can be understood from a local point of view 
based on first Hund's rule. We take C3& as examples. 
In Fig. 0(a), three spin-up (one non-degenerate-ui, one 
doubly degenerate) bands below Ep are derived mainly 
from Cr-3cZ orbitals. Five of six electrons on Cr (in these 



corresponding bands) occupy the majority spin states to 
yield Nj = 5. The sixth electron occupies the p-d hy- 
bridized, doubly degenerate but only half-filled m j-band 
yielding Nj = 1. Consequently, N=N|— N|=4, and hence 
/i = 4p; B . 

The present study predicts that linear chains of 
C„(TM) compounds (TM=Cr,Ti,Mn,Fe,Co) with spe- 
cific n can show half-metallic behavior with a diver- 
sity of spin-dependent electronic properties. Here, the 
type and number of atoms in the compound, as well 
as even-odd n disparity are critical variables available 
to engineer nanostructures with spin-dependent proper- 
ties. The electronic transport properties and the value 
of fi can be modified also by applied axial strain. Not 
only periodic structures, but also non-periodic combi- 
nations comprising HM-HM or HM-S (or M) quantum 
structures and superlattices can be envisaged to obtain 
desired device characteristics, such as spin-valve effect 
and spin-resonant tunnelling. Since linear carbon chains 
have been obtained also at the center of multi-wall car- 
bon nanotubes fig), C n Cr chains can, in principle, be 
produced inside a nanotube to protect the spintronic de- 
vice from the undesired external effects or oxidation. In 
fact, we obtained that strained C7C1" compound chain 
placed inside a (8,0) SWNT can be a HM. Of course, the 
properties revealed in this study correspond to idealized 
infinite chain structures, and are subject to modifications 
when the chain size becomes finite. However, for finite 
but long chains (for example a coil of C n TM around an 
insulating SWNT), the level spacings are still small to 
gain a band-like behavior. Also, localization of electronic 
states due to imperfections in ID may not lead to seri- 
ous difficulties when the localization length £ is larger 
than the length of the device. It is also noted that the 
properties of chains may depend on the type and detailed 
atomic structure of the electrodes. 

In conclusion, wc showed that half-metallic properties 
can be realized in linear chains of carbon-transition metal 
compounds presenting a number of exciting properties 
which can be of fundamental and technological interest 
for new generation devices. We believe that in view of 
recent progress made in synthesizing C atomic chains, 
present study will bring a new perspective in spintronics. 
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